Positions of Point-Nodes in Borocarbide Superconductor YNi 2 B 2 C 
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To determine the superconducting gap function of YNi2B2C, we calculate the local density of 
states (LDOS) around a single vortex core with the use of Eilenberger theory and the band structure 
calculated by local density approximation assuming various gap structures with point-nodes at 
different positions. We also calculate the angular-dependent heat capacity in the vortex state on 
the basis of the Doppler-Shift method. Comparing our results with the STM/STS experiment, 
the angular-dependent heat capacity and thermal conductivity, we propose the gap-structure of 
YM2B2C, which has the point-nodes and gap minima along (110). Our gap-structure is consistent 
with all results of angular-resolved experiments. 



PACS numbers: 74.20.Fg, 74.20.Rp,74.25. Jb 



I. INTRODUCTION 



The discovery of the nonmagnetic borocarbide super- 
conductor YM2B2C 1 has considerable attention because 
of the growing evidence for highly anisotropic supercon- 
ducting gap and high superconducting transition temper- 
ature 15.5 K. The boron isotope effect supports the clas- 
sification of this material as electron-phonon mediated 
superconductor^ At an early stage, from specific heat, 
thermal conductivity, Raman scattering and photocmis- 
sion spectroscopy experiments on YM2B2C, a highly 
anisotropic gap function was concluded^ In recent years, 
Maki et al. 5 - theoretically suggested that the gap sym- 
metry of this material is s+g wave and the gap function 
has zero points (point-nodes) in momentum space. Mo- 
tivated by this prediction, field-angle-dependent heat ca- 
pacity (FAD heat capacity and angular variation of the 
thermal transport (AV thermal transport) 7 - on YM2B2C 
have been measured and these results were considered to 
be consistent with this prediction. The FAD heat capac- 
ity and the AV thermal conductivity in H rotated within 
the ab plane show a fourfold oscillation with narrow cusps 
because of the presence of nodal quasiparticles subject to 
Doppler shifts. These results suggest that the gap func- 
tion has point-nodes along the a and b axes. However, in 
the analysis of both experiments, the isotropic Fermi sur- 
face (FS) was assumed. YM2B2C has highly anisotropic 
Fermi surfaces (FSs)*^^ 

The local density of states (LDOS) in the isolated vor- 
tex of YM2B2C was measured at 0.46K by the scanning 
tunneling microscopy and spectroscopy (STM/STS). 11 
The vortex core was found to be fourfold star-shaped 
in real space (see Fig. []Ja)-(c)). Near E — OmeV, the 
LDOS extends toward (100) (the a-axis) and has a peak 
at the center of the vortex core (Fig.[lja)). With increas- 
ing energy, the peak of the LDOS splits into four peaks 



toward (110) (Fig.[]Jb),(c)). The anisotropic spatial pat- 
tern of the LDOS around a vortex core is a consequence 
of the anisotropic pair-potential and FS. Many theoret- 
ical studies have been done since the LDOS of NbSe2 
can be observed by STM/STS.i2^1il5J&lLlS^2 i 2a We 
have calculated the LDOS pattern around a vortex core 
of YM2B2C, assuming that the gap symmetry is s+g- 
wave and the FS is isotropic^ Since this result explains 
the STM results only partly, obviously the calculations 
of the LDOS with a highly anisotropic FS are needed. 

The band structure for YM2B2C has been calculated, 
on the basis of the local density approximation (LDA). 
Lee et o/4 and SingbJi obtained three bands (17th, 18th 
and 19th) with mainly Ni-3d and Y-4d character which 
cut across the Fermi level by using the linearized muffin- 
tin orbital method and the general potential linearized 
augmented plane wave method, respectively. In recent 
years, band structure calculations were carried out by 
Yamauchi et al., 10 who are some of the authors, by us- 
ing a full potential LAPW(FLAPW) method. They have 
carried out LDA-based calculations with a certain modi- 
fication to describe well the experimental FSs which are 
obtained by de Haas van Alphen effect. In this modifi- 
cation, Y-d and Ni-d levels are shifted upward from the 
original LDA levels by O.llRy and 0.05 Ry, respectively. 
Such the modification have been successful in the FS cal- 
culations in LaB 6J 2i YbAl^ and LaRh 3 B 2 . ?3 This shift 
may originate mainly in the self- interaction and/or the 
non-local corrections to the LDA. They have shown that 
the density of states from 17th band, which is related to 
a multiple connected electron FS, has a sharp peak at 
Fermi energy. This peak comes from a van Hove singu- 
larity around (1/5, 1/5, 0) point in k-space as we will 
explain later in this paper. They have suggested that 
this singularity may lead large electron- phonon coupling 
locally, and give rise to anisotropic gap behavior in the 
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FIG. 1: (Color Online) Tunneling conductance images ob- 
served by Nishimori et al. at 0.46K in 0.07 T (72 x 72nm) 
for the bias voltage (a) OmeV, (b)-1.5meV, (c)-2.2meV. In 
this figure, we quoted from Fig. 5 (a), (c) and (e) in 'First 
Observation of the Fourfold- symmetric and Quantum Regime 
Vortex Core in YNi2B2C by Scanning Tunneling Microscopy 
and Spectroscopy' (H. Nishimori et al, Journal of the Physical 
Society of Japan 73, 3247 (2004)). The LDOS images calcu- 
lated for (d)e/A oc = 0, (e)0.15, (f)0.3, where 4v^£o X 4^2£ 
is shown. Here, Aoo denotes a pair-potential in the bulk re- 
gion averaged on the Fermi surface and £o denotes a coherent 
length. The smearing factor is S — 0.05Aoo 



superconducting state. 

The purpose of this paper is to determine the gap 
structure consistent with the STM/STS, the FAD heat 
capacity and the AV thermal transport experiments. We 
show the results of four calculations considering the band 
structure calculations. First, we calculate the LDOS pat- 
tern around a single vortex on the basis of the quasiclas- 
sical Eilenberger theory^ We assume various gap struc- 
tures, and finally obtain the gap structure which can ex- 
plain the STM experiment. Second, we calculate the den- 
sity of states in weak magnetic fields in the a-b plane on 
the basis of the Doppler-shift method^ We show that 
the FAD heat capacity and the AV thermal transport 
experiments are consistent with the angular dependence 
of the density of states from the gap structure consistent 
with the STM experiment in the presence of magnetic 
fields. Third, we show that the density of states from 
this gap structure is consistent with the STM experi- 
ment in the absence of magnetic fields. As results, we 



show the most suitable gap structure and positions of 
point-nodes. Finally, we show the LDOS patterns in the 
magnetic fields which tilts from the crystal c axis. With 
angular dependence of the LDOS, we can obtain more 
rich information about the gap structure. 



II. QUASICLASSICAL THEORY OF 
SUPERCONDUCTIVITY 



We consider the pair-potential; — 

A n (H,fe) = J dqA n (q,k)e i( l- R , (1) 



E 



(2) 

where k corresponds to the internal degree of freedom 
of the pairing state. Here, R is the center-of-mass coor- 
dinate of Cooper pairs, Vi l i tSltS2 (k, k') is the interaction 
between two electrons, 0g/ 2 ±fe' s * s ^ e aimi Wlation oper- 
ator for the quasiparticle states with spin s and momen- 
tum q/2±k' and we use units in which h = ks = 1- We 
assume the weak coupling interaction so that the pair- 
potential is not zero only near the Fermi surface. We con- 
sider pair-potential written as A(R, k) = iip(k)a y A(R) 
for singlet pairing. Here A(R) is a function of R. If 
the coherence length at the zero temperature £o is large 
compared to the Fermi wave length (~ 1/fcp), we can 
calculate the LDOS around a vortex core on the basis 
of the quasiclassical theory of superconductivity i 26 ' 27 i 28 
We consider the quasiclassical Green function g that has 
matrix elements in the Nambu (particle-hole) space as 



g(r,k,foj n ) = 



9 f 
-f 9 



(3) 



where uj n is the Matsubara frequency. The equation of 
motion for g called Eilenberger equation for singlet pair- 
ing is written as 



iv F (k) ■ Vg 



iw„ -A 
A* —fa} n 



■9 



= 0. 



(4) 



Here, vf is the Fermi velocity and [a, b] denotes the com- 
mutator ab — ba. The Green function g satisfies the nor- 
malization condition: g 2 = 1, where 1 is a 2 x 2 unit 
matrix. Considering clean superconductors in the type 
II limit we neglect the self-energy part of Green func- 
tion and the vector potential. 

The local density of states with the isotropic Fermi 
surface is given by 



Kr,e) = -K0) 



47T 



Re tr(<? r 



(5) 



Here, <7 R is the retarded Green function g R = g{\oj n — > 
e + iS) and v(0) denotes the density of states on Fermi 
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surface in the normal metallic state. The local density of 
states with the anisotropic Fermi surface is given by 



(6) 



Here, dS F is the Fermi-surface area element and v F is the 
modulus of Fermi velocity. 



III. RICCATI FORMALISM AND 
KRAMER-PESCH APPROXIMATION 

The Eilenberger equation can be simplified by intro- 
ducing a parametrization for the propagators that satisfy 
the normalization condition. Propagators are defined as 
P± = (1 =p g)/2, which were originally introduced in the 
studies of vortex dynamics i 29 ' 30 Using these propagators, 
we obtain the scalar equations expressed in Riccati for- 
malism as follows: 



v F ■ Va + 
v F ■ V6_ 



2u n a + 
- 2u>„b_ 



a + A*a A 
- 6_A6_ 



-A = 0, 
f A* = 0, 



where 



.9 

N 



-N 



(1 — a + 6_) 2ia + 
-2i&_ ' -(l-&_a+) 
l 



(1 + a+b 








(l + 6_a+)- 1 



(7) 

(8) 



(9) 
(10) 



Since these equations ([7]) and (jHJl contain V only through 
v F ■ V, they reduce to a one-dimensional problem on a 
straight line, the direction of which is given by that of 
the Fermi velocity v F . We consider a single vortex along 
the Z axis parallel to the crystal c axis. When we take 
the X axis on the a-b plane, the Y axis is determined 
automatically. We denote by Qm, &m and cm the unit 
vectors along X, Y, Z axis, respectively. We introduce 
the vector written as 



R = Xa 



M 



Yb 



M 



Zc 



(11) 



The origin R = is put on the vortex center. Because 
of a translational symmetry along the Z axis, the pair- 
potentials A(H, k), A*(R, k) do not depend on Z in the 
Riccati equations ([7]) and ©, and hence a+ and de- 
pend on R only through X and Y. As a result, v F ■ V 
in the Riccati equations can be replaced by 

d d 

v F ■ V -> v F ± ■ V = v x-q^ +Vy 9y' ( 12 ) 

with vx = v p • &m and vy = v F ■ 6m- Here, v F j_ is 
the vector perpendicular to the Z axis by projecting the 
Fermi velocity v F on the om^m plane. We introduce x, 
y, r by 



R — xv 

r = xv 



yu 



Zc 



M • 



(13) 
(14) 

(15) 



with 



cos sin 9 v \j am 



sm cos w v 



b 



M 



(16) 



Here, is the angle between om and the velocity v F ±. 
The resultant Riccati equations are then written as 



dx 
db. 



v FA _h(k F )^± + 2oj n a + + a 2 + A* - A = 0, (17) 



DpiWfcp)^ 2uj n b- - b 2 _A + A* = 0. (18) 



Here, v F ±_ and h(k F ) denote by 

/ dS F \v FA 



v F ± 



JdS ¥ 



(19) 



v FA _h{k F ) = ^Jv x (k F )+v Y {k F ). (20) 

An approximate expression for the quasiclassical Green 
function near a vortex core for low energy was analyt- 
ically obtained by Kramer and Pesch.— We call their 
approximation 'Kramer-Pesch approximation(KPA)' in 
present paper. Eschrig has also obtained the quasiclas- 
sical Green function near a vortex core in the first or- 
der of impact parameter and energy with use of the Ric- 
cati formalism* 2 ^ The method by Eschrig is equivalent 
to KPA. Using KPA, we can calculate the quasiclassical 
Green function around the vortex core in a low energy 
region (\aj n \ <S |Aoo|). Here, Aoo denotes a pair potential 
in the bulk region 



A, 



/ dS F limr^oo A(r, fe) 
JdS F 



(21) 



By expanding the equations {TTJ and (fT5|) in the first 
order of y and |a> n |,— we obtain the approximate solution 
as 

g{r, k F ; \u> n ) ~ - J Bl LJfll — , (22) 



2um 



am i 



where 



E w i(k F ) ~ 
F(x) = 



£ h{k F ) 

1 



v F ±h(k F ) 



AooA 2 ^), (23) 



dx'f(x'), (24) 



with £ — v F /(nA 00 ). Here, f(x) describes the spatial 
variation of pair-potential and /(0) = 0, lim x _>Q f( x ) = 
Aqo and A(fep) describes the variation of pair-potential 
in momentum space and MAX(A) = 1. The local density 
of states with the anisotropic Fermi surface is given by 



dS F v F ±\\\e 



~2\F(x) 



47t 2 £o|'Uf| 



-S{e-E ani ). (25) 



To consider the smearing effects, we approximate 
Eq. ([25} as 



v(r,e) 



dS F v F ±\X\Se 



-2XF(x) 



4kHo\vf\ (e - £ ani ) 2 + ^ ' 



(26) 



with the smearing factor S. Therefore, determining the 
distribution of the Fermi velocity h(k-p) from the band 
structure and the distribution of the pair-potential A(fep) 
from the gap structure, one can obtain the LDOS pattern 
around a vortex easily. 



IV. IN THE CASE OF YNi 2 B 2 C 

A. Local Density of States around a Vortex core 

We calculate the LDOS pattern around a vortex with 
using the band structure calculated by Yamauchi et al^ 
We consider the FS from 17th band and neglect the FSs 
from 18th and 19th bands since the values of the den- 
sity of states for 17th, 18th, 19th band at Fermi energy 
are 48.64, 7.88 and 0.38 states/Ry, respectively, so that 
the electrons of 17th band contribute dominantly to the 
superconductivity (See, Fig. 4 in Ref. fToh . The FS has 
two kinds of particular points on the planes k z = and 
k z = 0.5. We call the vectors at these points 'vector A' 
and 'vector B', respectively (See, Fig. [5]). The 'vector 
A' is the nesting vector and the 'vector B' is the vec- 
tor which connects the pair with the antiferromagnetic 
fluctuations. The Fermi velocities at these 'particular 
points' are in the (100) direction. Maki et al. suggested 
that the pair-potential at these particular points given 
above is strongly suppressed because of an instability in 
the particle-hole channel which strongly depresses the ef- 
fective potential for Cooper pairing^ We also investigate 
the positional relation between these particular points 
and point-nodes in momentum space to test Maki's sce- 
nario. 

We assume that point-nodes are around the points 
which have the strong antiferromagnetic fluctuations. We 
calculate the LDOS patterns assuming various gap struc- 
tures where positions of point-nodes are different. Com- 
paring with the STM experiments, we obtain the most 
suitable gap structure as shown in Fig. [51 In Fig. [21 
solid circles • and open circles o indicate the positions of 
point-node and the local minima, respectively. It should 
be noted that point-nodes are only at the 'vector A' and 
the gap functions do not have point-node at the 'vector 
B'. If the gap functions have point-node at the 'vector 
B', the LDOS patterns are not consistent with the STM 
experiments. 

In Figs. QJd)-(f), we show the calculated LDOS with 
this most suitable gap structure for several bias energies 
e. It is seen from Fig. [2(d) that the fourfold star centered 
at a vortex core extends toward (100) for e/Aoo = 0. As 
shown in Figs. QJe) an d (f), the peak of the LDOS splits 
into four peaks toward ( 110 ) with increasing energy. 
These LDOS patterns in Figs. QJd)-(f) coincide with the 
observation in Figs. QJa)-(c), respectively. If we assume 
the gap structure which does not have local minima at o, 
the calculated LDOS patterns are not consistent with the 
STM/STS experiments, since the peak of this calculated 
LDOS splits into four peaks toward ( 100 ) with increas- 




FIG. 2: The cross section including the F point of FS with 
the bet Brillouin zone for the 17th band. A solid circle • and 
an open circle o indicate the positions of point-node and the 
local minima, respectively. An open triangle A, box □ and 
diamond indicate the F point, Z point, X point, respec- 
tively. 



ing energy. If we assume the isotropic gap structure with 
this FS, the calculated LDOS patterns are almost circle. 
Our analytical theory in the previous paper™ shows that 
the LDOS around a vortex core consists of the contri- 
bution of the quasiparticles with momentum where the 
pair-potential is large on the FS. In other words, we can 
obtain the information of anti-node direction on FS from 
STM/STS experiments. From our calculations of LDOS, 
we can show that the gap amplitude at • and o are smaller 
than the half of the maximum gap. Therefore, we need 
comparisons with FAD heat capacity and AV thermal 
transport experiments, from which we can obtain the in- 
formation of gap-nodes. 



B. Angular Dependence of the Density of States in 
the Weak Magnetic Fields 

Comparing with the FAD heat capacity and the AV 
thermal conductivity experiments, we calculate the den- 
sity of states with the Doppler-shift method. 24 Without 
impurity scattering (superclean limit), the quasiclassical 
Green function is written as 



g(r,p;\w n ) 



iv F ■ v s 



y/(u> n +iv F -v s ) 2 + \A\ 2 



(27) 



Here, v s is the supercurrent velocity around a vortex. 
Therefore, the density of states at zero-energy v(e = 0) 
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FIG. 3: Angular dependence of specific heat C s oc v(e — 0, (f>) 
with use of the band-structure of YNi2B2C. cj> is the angle 
between magnetic field direction and the a axis. 
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FIG. 4: Quasiparticle density of states without magnetic 
fields, the smearing factor 8 = 0.05Aoo. 



is "written as 



v{e = 0, <j>) = ( 



Re- 



d{4>) 



(28) 



Here, d{cf)) = \v-p ■ v s ((j))\, the bracket means averaging 
over both Fermi surface and unit cell of vortex lattice 
and (j) is the angle between magnetic fields and the a axis 
on the a-b plane. Since the core states do not contribute 
to the specific heat in low magnetic fields, we neglect the 
spatial variation of the magnitude of the pair-potential 
and consider the spatial variation of the phase of the 
pair-potential around a vortex. Assuming d <C Aoo, only 
the quasiparticles around nodes in momentum space con- 
tribute to the density of states. This assumption is ap- 
propriate in weak magnetic fields. 

Figure [3] shows the angular dependence of specific heat 
C s oc v(e — 0, </>) calculated with use of the band- 
structure of YNI2B2C and the gap-structure consistent 
with the STM experiments (See Fig. [2]). Here, the mag- 
netic field is applied within the a-b plane at an angle 4> 
from the a axis. It should be noted that narrow cusps 
appear at <p — and <f> = n/2. This angular variation is 
consistent with the FAD heat capacity and the AV ther- 
mal transport experiments as shown in Fig. 2 in Ref. @ 
and Fig. 2 in Ref. 0, respectively. These angular vari- 
ations can be explained by the fact that quasiparticles 
(QPs) around the gap-nodes in momentum space con- 
tribute to the DOS and QPs traveling parallel to the 
vortex do not contribute to the DOS (Note that QP with 
vf parallel to the vortex line gives no contribution in 
Eq. (|28[) because v s is perpendicular to the vortex line 
and d — \vp • v s \ = for such a QP.). In other words, 
the narrow cusps appear at the direction of the Fermi ve- 
locity of QPs at the gap-nodes. In the case of YNi 2 B 2 C, 
the direction of the Fermi velocity at • is parallel to the 
( 100). 



C. Density of States without Magnetic Fields 



Comparing with the result of the STM/STS experi- 
ment in the absence of magnetic fields,— we calculate 
the density of states in the absence of magnetic field in 
the bulk region. In Fig. [4] we show the DOS with use 
of the band-structure and the gap-structure. Here, we 
assume the smearing parameter i5 = 0.05Aoo. Any sin- 
gularities like that at E = O.SAooin the case of s+g wave 
superconductivity 5 do not occur in our gap-structure. 
Our result is consistent with the density of states by the 
STM/STS experiment shown in Fig. 4(c) in RefUl 



D. Angular Dependence of the LDOS patterns 



One can obtain the three-dimensional information of 
the pairing symmetry from STM/STS experiments un- 
der the magnetic field with various directions. We have 
shown the LDOS patterns for the vortex parallel to the b 
axis about s+g wave superconductor in Fig. 14 in Ref. H2I 
We calculate the LDOS patterns on the a-b plane with 
various directions of magnetic fields, since the LDOS will 
be observed by the STM/STS experiments on the clean 
surface perpendicular to the crystal c axis. As shown in 
Fig. [5l these LDOS patterns are different qualitatively. 
The four peaks in the LDOS patterns for the vortex par- 
allel to the c axis disappear in the LDOS patterns for the 
vortex tilting from the c axis by angle 9 — n/4 within 
the a-c plane as shown in Fig.[5fa). Only the two peaks 
remain in the LDOS pattern for the vortex tilting from 
the a axis by angle <j) = 7r/4 and tilting from the c axis 
by angle 9 = ir/4 as shown in Fig. 03(b). 
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FIG. 5: (Color Online) Distributions of the local density of 
states for the vortex tilting from the a axis by angle 9 and the 
c axis by angle 0.(a) ((f), 9) = (0, tt/4), (b) (0, 0) = (tt/4, tt/4). 
The smearing factor is S = 0.05Aoo and the energy is 
e = 0.3Aoo. These LDOS patterns will be observed by the 
STM/STS experiments from the c axis. 

V. DISCUSSIONS 

We considered the 17th band and neglect the 18th and 
19th bands. By the angular resolved photo emission 
experiment^ the isotropic superconducting gap exists 
on the FS from 18th band and the FS from 19th band is 
not observed. Since the LDOS pattern from the isotropic 
superconducting gap is almost isotropic, the LDOS from 
18th band does not contribute to the four peak struc- 
ture of the STM/STS experiment. Comparing Fig. QJc) 
with Fig. QJf), the LDOS at the center by the theoretical 
calculation is quite smaller than that at a vortex center 
by the STM/STS experiments. Considering the FS from 
18th band, it seems that this LDOS at a vortex center 
becomes large. Therefore, the calculations with use of 
the FSs from 17th and 18th bands are needed for quan- 
titative comparison with STM/STS experiment. 

The results of STM/STS experiments around a vortex 
can not be explained by the 'two-gap model' since the 
strong anisotropy of the gap-structure makes the four 
peak structure of this experiment. Two isotropic gaps 
make two isotropic LDOS patterns even if the highly 
anisotropic FS is considered. The LDOS patterns with 
the model of point-nodes and that with the two-gap 
model are qualitatively different in the zero energy re- 
gion as shown in Fig. Efd). 

The previous analyses of the FAD heat capacity and 
the AV thermal transport experiments are insufficient 
since the isotropic FS is assumed, so that these analyses 
lead to the wrong conclusion that the gap function has 



point-nodes along the a axis and b axis. The informa- 
tion about nodes we obtain is the directions of the Fermi 
velocity at the nodes on the FS, since the Doppler-shift 
method is based on the fact that the energy has the term 
proportional to Vf • v s as shown in Eq. [28] Therefore, 
these experiments suggest that the Fermi velocity at the 
nodes on the FS is parallel to the a axis (the direction 
of nesting vectors). In the case of the anisotropic FS, 
the directions of the point-nodes are not found only from 
these experiments. In the only case of the isotropic FS, 
the directions of the point-nodes are parallel to the Fermi 
velocity at the nodes. 

We have shown that the point-nodes in momentum 
space are at the 'vector A'. These results suggest that the 
strong antiferromagnetic fluctuations suppress the super- 
conducting order parameter at these points^^ The rea- 
son why the point-nodes do not exist at the 'the vector 
B' is not clear yet. This is a future problem. 



VI. CONCLUSION 

In conclusion, we calculated the LDOS patterns around 
a vortex core, the angular dependence of the DOS in weak 
magnetic field, the DOS in zero fields in the bulk region 
and the LDOS patterns for the vortex tilting from the 
crystal c axis with use of the FS from 17th band obtained 
by the band-calculations, assuming various gap struc- 
tures where positions of point-nodes are different. Com- 
paring our theoretical calculations with the results of the 
STM/STS, the FAD heat capacity and the AV thermal 
transport experiments, we determined the gap-structure 
which is consistent with these experiments. The point- 
nodes are at 'the vector A', which is along (110) as shown 
in Fig. [5] We also showed that the previous analyses of 
the FAD heat capacity and the AV thermal transport 
experiments are insufficient since the isotropic FS was 
assumed. Considering the anisotropic FS, we showed the 
most suitable gap structure. This gap structure is con- 
sistent with the results of the FAD heat capacity, the AV 
thermal transport, the density of states without magnetic 
fields by the STM and the local density of states around a 
vortex core by the STM. We hope that our results will be 
tested by the STM/STS experiments with the rotation 
of the magnetic fields. 
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